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A realistic shell-model study is performed for neutron-deficient tin isotopes up to mass A = 108. 
All shell-model ingredients, namely two-body matrix elements, single-particle energies, and effective 
charges for electric quadrupole transition operators, have been calculated by way of the many-body 
perturbation theory, starting from a low-momentum interaction derived from the high-precision 
CD-Bonn free nucleon-nucleon potential. The focus has been put on the enhanced quadrupole 
collectivity of these nuclei, which is testified by the observed large B(E2;0~t — > 2^)s. Our results 
evidence the crucial role played by the Z = 50 cross-shell excitations that need to be taken into 
account explicitly to obtain a satisfactory theoretical description of light tin isotopes. We find also 
that a relevant contribution comes from the calculated neutron effective charges, whose magnitudes 
exceed the standard empirical values. An original double-step procedure has been introduced to 
reduce effectively the model space in order to overcome the computational problem. 

PACS numbers: 21.60.Cs, 23.20.Lv, 27.60.+j 


Light tin isotopes have been an interesting laboratory 
since the early 90s, when the experimental efforts toward 
the observation of the doubly-closed 100 Sn became a sort 
of search of the “Holy Grail”. This has led to a certain 
amount of data that have improved our understanding of 
the structure of neutron-deficient isotopes, providing also 
a challenging ground for shell-model calculations. As a 
matter of fact, the study of light tin isotopes opened the 
way to a new generation of realistic shell-model calcula¬ 
tions HH23, an approach that has then flourished in the 
last two decades. 

In the last few years, a renewed experimental inter¬ 
est has arised in studying these nuclei, especially with 
the help of intermediate-energy Coulomb-excitation ex¬ 
periments that are able to provide information on the 
electric quadrupole-excitation properties. 

In particular, in 2013 two papers reported about the 
measurement of the electric-quadrupole O) 1 ' -A 2^ tran¬ 
sition rate in 104 Sn. The first work was performed at 
GSI [2], and the measured value is B(E2-,0f —y 2+) = 
1000 ± 400 e 2 fm 4 . This value fits well with the predic¬ 
tions of realistic shell-model calculations [5], where em¬ 
pirical effective charges have been employed, both when 
the model space is made up by only neutron orbitals 
above the 100 Sn core and when also proton excitations 
coming from the proton Ogg /2 orbital are included con¬ 
sidering 90 Zr as a closed core. 

The second paper reported the results of an experiment 
carried out at the National Superconducting Cyclotron 
Laboratory (NSCL) at the Michigan University [6]. In 
this work the measured value of the 104 Sn B(E 2;0^ —► 
2+) is larger, 1800 ±370 e 2 fm 4 , and disagrees more than 
one sigma with the value of Ref. [3]. 

Very recently, another measurement of this transition 
probability has been performed at RIKEN [7], wherein 


the result has been obtained from absolute Coulomb ex¬ 
citation cross sections. The reported value of the 104 Sn 
B(E2;0f —2(j*”) is 1730±280 e 2 fm 4 , consistent with the 
result of Ref. [B] ■ 

The value obtained in the last two papers is quite large 
and is not reproduced by shell-model calculations, even 
when proton degrees of freedom are explicitly taken into 
account. In fact, various calculations HHB1ISP ,0 have 
been recently performed using 100 Sn and 90 Zr as closed 
cores, and they all predict B(E2) values too small for 
the neutron-deficient tin isotopes as compared to the ex¬ 
perimental ones. Moreover, no significant improvement 
is obtained by including also neutron excitations across 
the N = 50 shell closure, these excitations leading only 
to a slight increase of the B(E 2)s [4j. 

This background has been the main motivation to 
perform realistic shell-model calculations for neutron- 
deficient tin isotopes, using both the standard 100 Sn 
neutron-only model space and a larger one that includes 
Z = 50 cross-shell excitations with 88 Sr as an inert core. 
The main new elements of these calculations with respect 
to the previous ones are the inclusion of proton excita¬ 
tions from the lpi /2 orbital and the use of microscopic 
effective charges, as well as an original procedure to re¬ 
duce the large model space when considering 88 Sr as a 
closed core. 

Another motivation is to revisit a region that, as men¬ 
tioned before, has been for us the starting point for the in¬ 
vestigation of the reliability of realistic shell-model calcu¬ 
lations, going from the p-shell region up to nuclei around 
doubly-closed 208 Pb core jTTHTS] . 

In our shell-model calculation we start from the high- 
precision nucleon-nucleon potential CD-Bonn m, whose 
high-momentum repulsive components are smoothed out 
using the Ui ow -k approach [IS] so as to derive an effective 
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FIG. 1: (Color online) (a) Experimental l4H8l 1101 (red sym¬ 
bols) and calculated (black squares) excitation energies of the 
yrast J n = 2 + states and (b) B(E2\ 2 4 —> 0 4 ) transition rates 
for tin isotopes up to IV = 58, when using neutron effective 
charges (I) reported in Table I] Blue squares refer to calcula¬ 
tions with e^ mp = 0.5e. Shell model calculations have been 
performed using only neutron degrees of freedom (see text for 
details). 


hamiltonian H e g by way of the time-dependent pertur¬ 
bation theory mm- The chosen cutoff momentum is 
A = 2.6 fm _1 , and H c s has been calculated including 
diagrams up to third-order in V) 0 w-k- 

From the effective hamiltonian both single-particle 
(SP) energies and two-body matrix elements (TBME) 
of the residual interaction have been obtained |21j . and 
we have derived consistently the effective charges of the 
electric quadrupole operators at the same perturbative 
order. 

At first, we consider a model space spanned by the five 
neutron sdgh orbitals placed above doubly closed 100 Sn, 
so as to compare with previous realistic shell-model calcu¬ 
lations. The calculated excitation energies of the yrast 2 + 
states (E^+) and the B(E 2; 0^ —> 2+) transition rates for 
tin isotopes up to N = 58 are reported in Fig. [lland com¬ 
pared with the recent experimental data H® Tl0| . It can 
be seen that, while the experimental behavior of both ex¬ 
citation energies and B{E 2)s is reproduced, the observed 
quadrupole collectivity is underestimated by our calcula¬ 
tions, as evidenced by the fact that the predicted E*%s 

and B(E2-0+ -> 2+)s are larger and smaller, respec¬ 
tively, than the experimental ones. The B(E 2)s are too 
small notwithstanding the very large theoretical neutron 
effective charges e n , as reported in Table [j] It should be 
noted that the theory provides state-dependent effective 
charges, and the values relative to the low-lying (0g 7 / 2 ) 


orbital exceed the unity, being therefore quite different 
from the standard empirical one (0.5e). 

We stress that this is a result at variance with those ob¬ 
tained in other regions. Actually, in our previous works 
for nuclei above 40 Ca [52], 48 Ca [H], and 132 Sn [53] cores, 
starting from the same realistic potential, we have ob¬ 
tained effective proton and neutron charges close to the 
empirical values. This seems to indicate that for nuclei 
around 100 Sn relevant components of the real wavefunc- 
tion lie outside the chosen model space, which induces a 
large renormalization of the theoretical effective electric- 
quadrupole operator. 

For the sake of completeness, in Fig. [l] are also re¬ 
ported the results obtained with e n = 0.5e (blue squares), 
which do not differ too much from those reported in Ref. 
[6], where the effective interaction has been derived from 
the chiral N 3 LO [52] and NNLO [55] NN potentials. 

On the above grounds, we have considered a larger 
model space so as to take explicitly into account the 
Z = 50 cross-shell excitations of protons jumping from 
the lpi/ 2 , 0gg/2 orbitals into the sdgh ones. Within this 
large model space we have derived the effective hamilto¬ 
nian H { Jg, where the superscript indicates the number of 
proton (seven) and neutron (five) orbitals we have con¬ 
sidered. In Tables [T[n] are reported both theoretical neu¬ 
tron and proton effective charges, which are closer to the 
usual empirical values (e® mp = 0.5 4- 0.8e, e® mp = 1.5e). 
However, some of the state-dependent e n are close to 
unity, which may be considered an anomalous value with 
respect to the standard ones. 

As it could be expected, the enlargement of the model 
space provides wave functions closer to the real ones, even 
if the large e„s indicate that relevant components are still 
missing. 

The major difficulty with is that it cannot be diag¬ 
onalized for any tin isotope with up-to-date shell-model 
codes. One has then to overcome the computational 
problem finding some way to reduce the dimensions of 
the matrices to be diagonalized, and consequently make 
the shell-model calculation feasible. 

This problem has been also faced in Ref. [5] , where the 
shell-model calculations have been performed allowing up 
to 4 p — 4 h proton 90 Zr core excitations only. 


TABLE I: Neutron effective charges of the electric quadrupole 
operator E2 for the model space with 100 Sn (I) and 88 Sr (II) 


as cores. 


Tl'alaja Tlblbjb 

(alienII b) (I) 

(a||e n ||5> (II) 

0^7/2 0<?7/2 

1.20 

0.94 

0<?7/2 ld 5 / 2 

1.27 

0.96 

0(?7 / 2 ld 3 / 2 

1.19 

0.95 

W 5/2 1^5/2 

0.81 

0.94 

W 5/2 1^3/2 

0.83 

0.97 

W 5/2 2 Si / 2 

0.79 

0.79 

ld 3 /2 ld 3 / 2 

0.87 

0.96 

ld 3 /2 2S\/2 

0.85 

0.79 

Ohn/2 0/in/2 

0.78 

0.87 
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TABLE II: Proton effective charges of the electric quadrupole 
operator E2. 


Tlalaja 

riblbjb 

{a\\e p \\b) 

O 

to 

CO 

to 

059/2 

1.62 

059/2 

057/2 

1.67 

059/2 

14/2 

1.60 

0 

to 

-sj 

to 

057/2 

1.73 

057/2 

1^5/2 

1.74 

057/2 

1^3/2 

1.76 

1^5/2 

1^5/2 

1.73 

1^5/2 

1^3/2 

1.72 

W5/2 

2S1/2 

1.76 

1^3/2 

1^3/2 

1.74 

1^3/2 

2S1/2 

1.76 

0^11/2 

0/lll/2 

1.72 


In the present work, we have resorted for the first time 
to an approach which, by way of a unitary transformation 
of leads to a new effective hamiltonian defined in a 
truncated model space. The choice of the truncation of 
the model space is driven by the behavior, as a function 
of Z and N, of the proton and neutron effective single¬ 
particle energies (ESPE) of the original hamiltonian 
so as to find out what are the most relevant degrees of 
freedom to describe the physics of light tin isotopes. To 
this end, we report in Figs. [2] and [3] the evolution of both 
proton and neutron ESPE as a function of Z. From the 
inspection of Fig. [2j it can be observed that an almost 
constant energy gap provides a separation between the 
subspace spanned by the lpi/ 2 , 0 . 59 / 2 , 1 ^ 5 / 2 , 057/2 proton 
orbitals and that spanned by the 2si/ 2 , ld 3 / 2 , 0 h 33 / 2 ones. 
This leads to the conclusion that a reasonable truncation 
is to consider only the lowest four orbitals, as proton 
model space. 

On the neutron side, Fig. [3] evidences that the filling 
of the proton Ogg / 2 orbital induces a relevant energy gap 
at Z = 50 between the 1 ^ 5 / 2 , 057/2 subspace and that 
spanned by the 2si/ 2 , ld 3 / 2 , 0/in / 2 orbitals. This gap, 
around 2.4 MeV, traces back to the tensor component of 
the proton-neutron interaction that is mainly responsible 
for the shell evolution [26]. Our calculated monopole 
component of the proton-neutron 0 g 9 / 2 , 057/ 2 interaction 
is -0.479 MeV. 

We, therefore, have deemed it reasonable that a neu¬ 
tron model space spanned only by the 1 ^ 5 / 2 , 057/2 or¬ 
bitals may provide the relevant features of the physics of 
light tin isotopes. 

Moreover, if we consider the evolution of the neutron 
ESPE as a function of N (see Fig. [4]), it can be observed 
that 2 S 1 / 2 , 1 ^ 3 / 2 , 0 ^ 11/2 orbitals start to play a more rel¬ 
evant role from 108 Sn on. Actually, the reduction of the 
energy gap between these orbitals and the 1 ^ 5 / 2 , 057/2 
ones, and the progressive filling of the latter, may reduce 
the effectiveness of the truncated model space to describe 
the physics of tin isotopes above N = 56. 



Z 

FIG. 2: Calculated proton effective single-particle energies of 
Hjg as a function of the atomic number Z. 
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FIG. 3: Calculated neutron effective single-particle energies 
of H as a function of Z. 


On these grounds, we have derived a new effec¬ 
tive hamiltonian H[fg, defined within a model space 
spanned only by the lpi/ 2 , 059 / 2 , lds/ 2 , O 57/2 proton and 
1 ^ 5 / 2 , 057/2 neutron orbitals, by way of a unitary trans¬ 
formation of Hjg (see, for example, Refs. mm)- 
It should be pointed out that we have applied this uni¬ 
tary transformation to the two valence-nucleon systems 
only, so that the energy spectra of 90 Zr, 90 Sr, and 90 Y 
are exactly the same when diagonalizing H and H^g. 
In order to obtain the same outcome for the eigenvalues 
of 102 Sn, besides one should include effective many- 
body forces, that can only be obtained by diagonalizing 
Hjg for this nucleus. As already pointed out, this is un¬ 
feasible, so we have taken into account only TBME of 
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FIG. 4: Effective single-particle energies of tin isotopes as a 
function of N calculated with Hjg. 



FIG. 5: Experimental |29' and theoretical spectra of 1(,2 Sn 
calculated with H^g and H^g (see text for details). 


It is obvious that the larger is the chosen subspace the 
smaller is the role of these effective many-body compo¬ 
nents. First of all, it may be worth verifying the relia¬ 
bility of our truncation scheme for tin isotopes. To this 
end we have derived another hamiltonian, starting from 
ifjg, defined in the full sdgh neutron model space but 
with the same proton model space of H'^'g. This is the 
largest model space in which we can manage to diago¬ 
nalize the shell-model hamiltonian of 102 Sn. We dub this 
effective hamiltonian H^g, and in Fig. p4we compare the 
low-energy spectra of 102 Sn obtained by means of H^g 
and H^g with the observed one [2$ ]. 


It can be noted that is able to reproduce quite 
well the experimental spectrum, and that the spectrum 
calculated with is in a good agreement with results 
obtained with H^g. It is very relevant, for the subject 
of our study, to point out that the B(E 2; Oj 1 " —>• 2+) cal¬ 
culated with H^g is 1065 e 2 fm , very close to the value 
of 1135 e * 2 * fm 4 obtained with H^g- This supports the 
adequacy of our truncation scheme when using H^g and 
increasing the number of valence neutrons. This is an im¬ 
portant result since our calculations, performed by way 
of the Oslo shell-model code [30], cannot be extended to 
heavier tin isotopes using H^g. 
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FIG. 6: (Color online) Same as in Fig. [I] but with shell-model 
results obtained with H*g. Blue squares refer to calculations 
with e e n mp = 0.5, e e p mp = 1.5. 


In Fig. jfij we report the experimental E^+ and 

B(E2-0f —► 2f) (red symbols) for the tin isotopes up 
to N = 58, and compare them with the results obtained 
using H^g (black squares). We have also included our 
prediction for the closed shell 100 Sn, where only proton 
degrees of freedom are taken into account. As in Fig. 
[T] we have included the calculated B(E2 )s when using 
e° mp = 0.5e, e; mp = 1.5e. 

We see that the shell-model calculations with H^g 
are able to reproduce quite well the experimental 

2 i 

and, when employing the theoretical effective charges, 
the B(E2;0f — 1 - 2^) up to A = 106 and, consequently, 
the onset of collectivity from 102 Sn on, driven by the 
Z = 50 cross-shell excitations. As a matter of fact, the 
wavefunctions evidence a depletion of the proton 0 <?g /2 
orbital from 102 Sn to 106 Sn, as testified by the occupa¬ 
tion numbers reported in Table |III| 

It is worth to point out that the agreement with ex- 
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periment for 108 Sn deteriorates with respect to lighter 
isotopes, owing to the the fact that, as already men¬ 
tioned, the influence of neutron ld 3 / 2 , 2si / 2 ,0/in / 2 or¬ 
bitals starts to play a non-negligible role. 


TABLE III: Occupation numbers of proton 

lPi/ 2 , 039 / 2 , 0 ^ 7 / 2 , ld 5/2 of 102 ~ 108 Sn J = 0^,2^ state, 

calculated with model space III (see text for details). 


Orbital 

1U2 Sn 

104 Sn 

iUo Sn 

iU8 Sn 

+ 

o 

II 

vrlpi /2 

1.97 

1.98 

1.98 

1.98 

7r0g 9 / 2 

9.55 

9.38 

9.36 

9.38 

71-0^7/2 

0.26 

0.28 

0.27 

0.25 

7lTd 5 /2 

0.22 

0.37 

0.40 

0.39 

to 

+ 

7rlPi/2 

1.98 

1.98 

1.98 

1.98 

7r0g 9 /2 

9.50 

9.23 

9.21 

9.27 

71-0^7/2 

0.25 

0.28 

0.27 

0.24 

7rld 5 /2 

0.28 

0.51 

0.54 

0.51 


In summary, we have performed a shell-model study of 
light tin isotopes starting from a realistic NN potential, 
where all results have been obtained without resorting to 
any empirical parameter. The main features of present 


work may be itemized as follows: 

• We have confirmed the crucial role of the Z = 50 
cross-shell excitations to obtain a satisfactory de¬ 
scription of tin isotopes. This implies the use of a 
large shell-model space including both proton and 
neutron orbitals. 

• We have followed an original double-step approach 
to reduce the computational complexity of the 
shell-model problem. This is based on the study 
of the ESPE of the large-scale hamiltonian, so as 
to identify the most relevant degrees of freedom 
to be taken into account in the construction of a 
truncated shell-model hamiltonian. To this end, a 
unitary transformation is employed. 

• We have highlighted the role of theoretical effective 
charges in reproducing the ciuadrupole collectivity 
of the B{E2)s. 

• We have presented some predictions for 100 >i 02 g n 
spectroscopic properties, which may provide guid¬ 
ance for future experiments. 
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